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The compatibility of the QCD sum rules and effective hadronic models predictions 
, are examined. For this purpose we have considered the results for the nucleon self- 

■ energy in a dense hadronic environment provided by two independent QCD sum- 
^\ . rules calculations. They are immersed in a theory of hadronic fields giving rise to 

non- linear interactions, whose vertices are parameterized in different ways. Although 
all of them reproduce the self-energy used as input, very different descriptions of 

■ nuclear observables are obtained. Only under very definite circumstances we have 
found an acceptable agreement with the nuclear matter properties. To achieve this, 
phenomenological parameters are not required at all. 
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\ Currently, there is a noticeable interest in interpreting and formalizing the low energy 

manifestations of the fundamental theory of the strong interactions. However, the mathe- 
matical complexity that Quantum Chromodynamics (QCD) exhibits in this regime makes 
impossible to reach this aim. Several approaches were proposed to circumvent this dif- 
ficulty, such as lattice simulations, QCD sum-rules, and effective models, for instance 
Nambu- Jona Lasinio, bag-like models, and chiral perturbation theory. Each of these treat- 
ments emphasizes some aspects of QCD, considered as the most relevant in the hadronic 
phase. Contact with low energy physics is ensured by fixing some model parameters, 
such as hadronic masses and decay constants to its experimental values. It is generally 
believed that for high density and/or temperatures, the symmetries of the fundamental 
theory would be recovered, therefore effective models should reproduce this feature to 
some degree. 

The low to medium density states of nuclear matter have been extensively studied by 
using nuclear potentials. Most of the nuclear structures, reactions, and dispersion phe- 
nomenology have been successfully described in this way. However, these formulations lack 
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of Lorentz covariance and they violate causality when extrapolated to higher densities. 
The requirements of covariance and causality can be fulfilled properly within the hadronic 
field theory (HFT), whose prototype the Quantum Hadrodynamics [1] usually takes dense 
nuclear matter as its reference state. This practice contrasts with certain nuclear potential 
treatments that refer to in- vacuum scattering lengths, ascribing this feature to a more 
fundamental instance of nuclear potentials as compared with HFT. Inspired in this view- 
point, several investigations have been made in order to make compatible both schemes, 
due to the simplicity and accuracy of the HFT treatment. The in-medium self-energy of 
the nucleon plays a key role in these attempts, as it is evaluated in both approaches and 
then equated to extract medium dependent vertices [2,3,4,5]. 

The idea behind this procedure is that when a piece of information provided by another 
field of research is coherently inserted into the framework of the HFT, a simpler and more 
intuitive description is obtained instead of stating first principles interactions and per- 
forming involved calculations. This procedure is specially suited for high density systems 
since the details of the interactions are faded out in favor of a statistical average. A signi- 
ficative exemplification is given by the Brown- Rho scaling of hadronic masses. Taking into 
account the chiral and scale symmetries of QCD solely, an approximate scaling law for the 
in-medium hadronic masses was derived in [6]. This hypothesis was applied to describe 
heavy ion collision, reaching an excellent agreement with the experimental results for the 
low mass dilepton production rate [7]. 

On the contrary, going from low to higher energies, references [2,3,4,5] take nuclear poten- 
tials as the starting point and results are extrapolated to the medium-high dense nuclear 
matter and atomic nuclei. In ref. [4] it was stressed that a coherent model requires that 
the medium dependence of the above mentioned vertices must be introduced in terms of 
hadronic fields. Otherwise, the rearrangement self-energies are missed and thermodynam- 
ical consistency is not satisfied. 

In this work, we try to match QCD inspired results with the HFT formalism, by intro- 
ducing field dependent strength couplings as proposed in [4]. For this purpose, we have 
chosen as inputs the nucleon self-energies evaluated within the QCD sum-rules approach 
[8,9,10,11,12]. We have compared different field parameterizations of the vertices, and ex- 
amined the ability of our results to fit the nuclear matter phenomenology. 
QCD sum-rules is a powerful procedure to reveal the QCD foundations of certain hadronic 
properties. This method relies on the evaluation of correlation functions in terms of quarks 
and gluons degrees, then applying the operator product expansion it is possible to express 
them as combinations of perturbative contributions and condensates (non-perturbative) . 
In the last step, results are connected with the hadronic counterpart, from which infor- 
mation on the spectral function can be inferred. The method was developed to study 
meson [8] as well as baryon [9] properties in vacuum, it was subsequently generalized to 
study finite density systems [10,11,12]. We have used the in-medium nucleon self-energies 
provided by references [11,12], more precisely the simplified expressions given there as 
functions of the baryonic density. Both investigations include non-perturbative contribu- 
tions representing one meson exchange, the structure of the nucleon, and multiple meson 
exchange. Nevertheless, they differ in the criteria used to simplify the final formulae. The 
results of ref. [12] keep all the contributions averaged over the Borel mass parameter; 
whereas, in [11] only the leading term in the operator product expansion is retained, and 
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the Borel mass is fixed in order that the effective nucleon mass goes to its experimental 
value in vacuum. In the latter linear dependence on the density is obtained for 

both the scalar and vector nucleon self- energy. 

As it was above mentioned, the Density Dependent Hadron Field Theory (DDHFT) [4,5] 
attempts to insert the nucleon self-energy supplied by Dirac-Brueckner calculations with 
one boson exchange potentials, into a hadronic field description of the nuclear matter 
phenomenology. In general terms, we have adopted the procedure described in ref. [4] to 
merge dynamically the external input into a field model, and we will indicate where we 
depart from it. 

We start with a lagrangian for nucleons (V') and scalar and vector mesons {a and uj 
respectively), resembling the QHD-I model of Walecka [1] for isospin symmetric nuclear 
matter, 

jC^'4^{i^-M + Gsa-G,i/>)iP + ^{d^ad'^a - my) - ^F^.F'^'^ + ^m>^u;^ (1) 

where F^^, = d^uj^ — d^uj^ has been used, and the meson masses have been fixed at the 
empirical values = 783 MeV, = 550 MeV. The main difference is that Gg and G^ 
are not coupling constants as in the Walecka model, but they are field functionals. Since 
these vertices are Lorentz scalars, we expect that the dependence on the fermion fields 
appears through the bihnear combinations si(x) = ■ijj{x)ip{x), S2{x) = ^Jn^J,{x)n^^{x), with 
^^(3^) = V'(^)7mV^(^) ^he baryonic current. Ref. [4] only developed the dependence on Si 
or S2 ; in this paper, we also consider the eventual dependence on the meson fields but 
neglect derivatives of any of the fields, namely Gs^v = Gs^y{si, S2, S3, S4), with S3 = ^cu^cu^, 
and S4 = a. The field equations are obtained in the usual way, by minimizing the action 
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duji, \ ds 



dG. , dG, ; 



(2) 
(3) 
(4) 



The two last terms within brackets in Eq. (2) evaluated in the mean field approximation 
can be regarded as the scalar and vector components of the nucleon self-energies. In order 
to make the problem solvable, only a set of variables can be taken into account in the 
evaluation of the functionals Gg and G^. We have considered the following cases: a) both 
vertices depend only on S2, b) Gg depends on si, G^ depends on S2, and c) Gg depends 
on (T, Gy depends on S3. The case (a) has been developed in great detail in ref. [4], the 
same as case (b) they produce meson equations resembling those of the Walecka model, 
and nucleon self-energies with the rearrangement contributions added. In the instance (c) 
the self-energies have the same structure as in QHD-I, but the source term of the meson 
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equations are modified. 



Tlic mean field approximation (MFA) is suited to describe static homogeneous nuclear 
matter. Within this scheme, meson fields are replaced by their uniform mean values and 
bilinear combinations of spinors are replaced by their expectation values. Therefore within 
MFA, Eqs. (2-4) reduce to 



rriga — < ip Ggip > + < ip-;: — a > 



2~v 



da 

< -0 -0 > - < V'^r^ — V' > + < V'^r^ — > 



(5) 
(6) 

(7) 



where tildes over the meson symbols stand for their mean field values, and 



QQ QQ 

j:,^<Gsa> + <ij^-^a-ij^x'^p-Qfi^^>, (8) 

1 - dG - dG 

= < G^ou^ >+-< {-ijij^a + ijjxi'^'^^)'ip%'ip > ■ (9) 

S2 dS2 OS2 

These results are simplified within the assumptions (a-c) above, for the sake of complete- 
ness we consider separately each of these cases: 



a) Hs^GgCr, E^ = G^u^ (dG^/dn) u'^A^^ - {dGjdn) ctA^, a = Gsnjml, 
a;^ = G^n^/m^,with A^^ =< ■4)^ui^i)^^j,'il^ > /n, and A^ =< 'ijjijj'ijj^^ijj > /n, 

b) Es = 5-(G's + Us dGs/dus), E^ = GyU/j, + (dGy/dn) cL'^'A^^, 

a, and cD^have the same expressions as in (a), 

c) Eg^Gsd-, E^ = GyUJf,, mla = GgUs + ans{dGs/da), 

mluj^ = GyUij, {dGy/du}) n^Cj''uj^j,/ sz, 

we have used Gg^v —< Gs,v >, —< "ipip >, =< ■07//'0 >, n — Uq is the baryonic 
density, and we have taken into account the fact that S2 = n, as it is clear when it is eval- 
uated in the reference system of rest matter. Finally, we have introduced the simplifying 
assumption that Gs,v and its derivatives can be extracted from the expectation values, 
even when they are field functionals. Regarding to the expectation values A^^, and A^, 
they can be evaluated by using the Wick theorem and the propagator for the nucleons 
in the medium, giving rise to direct and exchange terms. In our calculations we have 
neglected the Fock terms, in which case we have obtained A^,^ = rii/n^/n, A^ = n^n^/n. 

In previous investigations [2,3,4], the procedure to define the functionals Gg and Gy con- 
sisted of including the nucleon self-energies as input into the equations 

^fP = Ggns/ml E7 = G,n^/m^, (10) 
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where S*"'^, E^"^ are the self-energies obtained with one boson exchange potentials in 
the Dirac-Brueckner-Hartree-Fock approach for symmetric nuclear matter. The Eqs. (10) 
coincide with Eqs. (8)- (9) if the rearrangement terms were omitted, in spite of the impor- 
tance of these contributions for satisfying the thermodynamic consistency, as stressed in 
ref. [4]. At this point we deviate from the procedure of the DDHFT in two aspects; firstly, 
we adopt as input the nucleon self-energies 2*^*^°^^ grounded on the fundamental theory 
of strong interactions. Secondly, we insert this input into the full expressions (8) and (9), 
in fact we regard from the very beginning the field dependence of Gs^v Therefore, the 
vertices are the solution of a differential equation which must be solved together with the 
self-consistent condition for a in cases a) and b), or a pair of differential equations for a 
and Cj in the case c) . The explicit form of these solutions are shown below, for each of the 
instances considered 




where we used that the spatial components of tetravectors become null for isotropic mat- 
ter, and we considered E^^^^'^ as being parameterized in terms of only the baryonic density 
n. The expressions for Gs,v become singular at zero density, nevertheless they have finite 
limits: mldll^'^'°^'^{0) / dn, Gl mldE'^^°^^{0) / dn, for all the cases under consid- 

eration. 

Up to this point, we have constructed a hadronic effective interaction capable of reproduce 
the nucleon self-energies as they were evaluated in the QCD sum-rules; therefore, our 
findings are subjected to the same limitations in applicability. 

In ref. [14], a similar purpose was attained by adopting the results of [11] and simulating 
the hadronic interaction by means of four fermion couplings. The vertices are taken as 
functions of only they were determined in the same way as in the DDHFT, namely, 
using Eq. (10). Nevertheless, the authors consider that the effect of chiral fluctuations 
is predominant in the many body aspects of nuclear matter, and the variation of the 
self-energies according to QCD sum-rules calculations has a stabilization effect at high 
densities. 

In the next step we have evaluated the energy-momentum tensor T'^'^ by the canonical 
procedure, and finally, the energy per unit volume in the MFA is obtained by taking the 
in-medium expectation value of T^. This gives 

PF ,3 

Emfa^ I jA-Jp' + M*^ + n,(Ef °- - G,a) + -{mla' + mlu:'), (11) 
^ ' 
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where we have introduced the Fermi momentum pp, related to the baryonic density by 
n = 2p'p/{37i^), and the effective nucleon mass M* = M — S^*^"^^. It must be noticed 
that the sign of the scalar self-energy has been changed as regards the one used in refs. 
[11,12], in order to adopt the practices of QHD models. 

Using the thermodynamical relation P — /in — Em fa we have obtained the pressure of 
the nuclear matter, with the chemical potential given hj ^ — Ep + for cases a) and 
b), and ji = Ep + ndTiy/dn for the remaining one, being Ep — •y'p|r+M*2. 
A feature of isospin symmetric nuclear matter is that the energy per particle has a mini- 
mum, giving rise to bound states at zero temperature. This is the least requirement that 
a model of nuclear matter should satisfy. A measure of this property is the binding energy 
Eb — EMFA/n — M, which is well determined to have a minimum value Eb ^ —16 MeV 
at the normal density Us = 0.15/m~^. We have examined the binding energy in the three 
cases studied for the functional dependence of Gg^v, which are displayed in Fig. 1 in terms 
of the baryonic density. As we can see, none of these instances is able to fit the saturation 
properties of nuclear matter. Our results with the input taken from ref. [11] does not 
distinguish among the different parameterizations of the vertices; on the other hand, the 
case c) evaluated with the self-energies given by ref. [12] (case cD in the following) clearly 
stands out from the other ones. The non-monotonic behavior of the curve cD is more 
compatible with the phenomenology, more precisely the binding energy has a minimum 
Eb — —4 MeV at the density n — 0.8 n^, not far from the empirical saturation point. 
It is apparent that the approaches a) and b) yield undistinguishable results for Eb- On 
the other hand, the case c) produces higher meson mean values, enhancing the energy 
density. This fact reverses the decreasing trend of Eb obtained for the ) and b) in 

the right hand side of this figure, the net behavior is non-monotonous as shown there. On 
the other hand, the same effect applied to the increasing curves (a-b) on the left panel, 
reinforces the growth. Although, the small increment is not appreciable within the scale 
of this figure. 

It has been stated that in order to obtain the slight binding energy and simultaneously 
a high nuclear spin-orbit potential, the scalar and vector components of the nucleon self- 
energy must have nearly the same value, of roughly three hundred MeV for densities about 
Us [13]. This requisite is satisfied by both refs. [11] and [12]; however, as can be deduced 
from the results a) and b) it is not sufficient to determine conclusively the many body 
aspects of the system. 

As the QCD sum-rules calculations are subjected to uncertainties of diverse origins, we 
have introduced two independent scaling factors Xg and Xy for each of the nucleon self- 
energies used as inputs, and we have explored the range | Xs,v — 1 |< 0.25. We have as- 
sumed these factors independent of the density, consequently their derivatives are scaled 
by the same magnitude. We have found that within this range of variation, only the case 
(cD) is able to reproduce simultaneously the binding energy and the density of satu- 
ration with a slight scaling Xg ~ 1.19, Xy ~ 1.12. Aiming to a quahtative description 
for all the densities considered, we have introduced a uniform scaling of the self-energies 
E^°°®'^(n) Xs^vT,'^^°^^{n) within the scheme c). 

As the next step, we have examined the behavior of the Landau parameters Fq and Fi 
for symmetric nuclear matter. They can be evaluated as the Legendre projections of the 
second derivative of Emfa with respect to the baryonic density, for a deduction see [15]. 
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These parameters are very useful in regard to collective phenomena of the dense nuclear 
environment, as for instance, the giant monopolar and quadrupolar modes, or the sound 
velocity [16]. The nuclear compressibility can be evaluated as = 3pp{l + Fq)/Ep, its 
value at the equilibrium density is constrained experimentally to range between 230 MeV 
< K < 270 MeV. Within the approach (cD) we have obtained K = 245 MeV, and re- 
lated to it is the giant monopole mode energy Em — cuq/A^/^. In our calculations we 
have obtained uq = 108 MeV, whereas, the experience yields c<Jq^^ = 80 MeV. On the 
other hand the scaled calculations lead to the excessively high values K = 660 MeV, and 
cuq — 180 MeV. Another interesting property is the excitation energy of the quadrupole 
state, given by the relation Eq — 002/ A^^^ with 002 — PfV^/ (1-2/iy^l + Ei/3), yielding in 
our calculations U2 — 78 MeV, and 102 — 100 MeV in the non-scaled and scaled cases, 
respectively. These values must be compared with the empirical result 002^^ = 63 MeV. 
On the other hand, the numerical values obtained for c^o and 002 without scaling, can be 
well accommodated among the theoretical results for the isoscalar giant resonances in the 
non-hnear models of HFT [17]. Finally, the sound velocity Vg = Pp{l + Fo)/{3hEf) takes 
on the value Vg/c — 0.17 and Vg/c — 0.28 at the saturation density for the non-scaled 
and scaled instances, respectively. The magnitude of Vg increases with density quickly, 
becoming acausal before n ~ 2ns is reached. This feature puts severe limits to the range 
of applicability of the scheme. 

In Fig. 2 we show the vertices Gg^y for the (cD) approach, in terms of their respective 
variables. They are well fitted by the expressions Gg = a/{l + bss) and Gy — A + B S4, 
with a = 7, 6 = -0.59fm-\ A = 7.83, and B = 8.39/m-^ 

To sum up, in the present work we have examined the possibility of matching QCD sum- 
rules predictions and hadronic field theories. A successful combination would result in a 
profitable feedback, by giving foundations to hadronic models and corroborating simul- 
taneously the QCD sum-rules methods and findings. In particular we have used results 
for the micleon self-energy in a dense hadronic environment derived in two independent 
treatments of QCD sum-rules, with diverse degree of approximation. As a first step, we 
have studied different formal schemes of merging these inputs into hadronic field models 
and we have checked their abihty to fit the basic nuclear matter phenomenology. We have 
introduced interaction vertices Gg and Gy, which are functionals of the hadronic fields, 
and we have determined them by requiring the reproduction of the self-energies used as in- 
put. It must be emphasized that in this procedure we have included the full expression for 
the self-energies, having obtained a set of differential equations for the definition of these 
vertices, instead of algebraic ones,. We have found that only a parameterization of the 
vertices in terms of mesonic degrees is compatible with the basic features of dense nuclear 
matter. In addition we have presented a fitting in terms of rational functions, which make 
evident the non-linear quality of the meson-nucleon coupling. If we allow an uncertainty 
margin of 20% in the sum-rules calculations, this scheme is able to fit precisely the sat- 
uration properties of the isospin symmetric nuclear matter. Another physical quantities, 
which are related to collective nuclear phenomena, agree with the order of magnitude of 
the experimental values, although they are not admissible as a quantitative description. 
It must be stressed that no adjustable parameters are introduced, only the QCD sum-rules 
self-energies are required as input. The parameters Xs,Xy were introduced only to study 
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the versatility of the results. 



We conclude that compatibility of QCD sum-rules and hadronic field models is feasible. 
However, a detailed qualitative description will require further improvements in QCD 
sum-rules calculations, in particular the high density QCD phenomenology should be 
properly included. Therefore, additional refinements and higher order calculations will be 
useful and deserve to be investigated. Improved results could be inserted into hadronic 
models to test another significative manifestations of the hadronic phase such as finite 
nuclei and the eventual deconfinement phase transition. 
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Fig. 1. The binding energy as a function of the baryonic density for each of the paramctcrizations 
of the interaction vertices, cases a-c. By using as input the self-energy given by ref. [11] all the 
curves coincide, in the case of including the nucleon self-energy of ref. [12] the case c) distinguishes 
clearly from the other ones. The line convention is indicated in the figure. 
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Fig. 2. The interaction vertices of the hadronic field model as functions of the meson mean field 
values. The model reproduces the nucleon self-energy given by ref. [12]. 
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